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Abstract 

We consider the existence of localized modes corresponding to eigenvalues of the periodic 
Schrodinger operator —d^ + V{x) with an interface. The interface is modeled by a jump ei- 
ther in the value or the derivative of V{x) and, in general, does not correspond to a localized 
perturbation of the perfectly periodic operator. The periodic potentials on each side of the 
interface can, moreover, be different. As we show, eigenvalues can only occur in spectral gaps. 
We pose the eigenvalue problem as a gluing problem for the fundamental solutions (Bloch 
functions) of the second order ODEs on each side of the interface. The problem is thus reduced 

to finding matchings of the ratio functions R± = ^^^^^ , where '0± are those Bloch functions 
that decay on the respective half-lines. These ratio functions are analyzed with the help of the 
Priifer transformation. The limit values of R± at band edges depend on the ordering of Dirichlet 
and Neumann eigenvalues at gap edges. We show that the ordering can be determined in the 
first two gaps via variational analysis for potentials satisfying certain monotonicity conditions. 
Numerical computations of interface eigenvalues are presented to corroborate the analysis. 

1 Introduction 

Localization for perturbed periodic Schrodinger operators L = —A -|- Vo{x) + V{x), where Vo{x) 
is periodic in x G M",n G N, is a classical problem traditionally treated by spectral theory. Most 
commonly it is studied for perturbations V{x) that are either compactly supported, see, e.g., Deift 
& Hempel [6], Alama et al. [1 , and Borisov & Gadyl'shin [3^ or fast decaying, e.g. V £ L"'/^(]R"), 
cf. Zeludev [23] and Alama et al. |T] . Both of these scenarios can lead to eigenvalues of L and thus 
to localization. Potentials V describing random perturbation also yield eigenvalues due to Anderson 
localization, studied, for example, by Kirsch et al. [12j and Veselic [21|. We investigate localization 
in the one-dimensional case n = 1 due to the presence of deterministic interfaces which cannot 
be represented as localized perturbations of —d'^ + Vo{x). Such an interface arises, for instance, 
when V{x) is periodic on one side of the interface and vanishes on the other side (we assume 
commensurability of the periods of V and Vq to preserve periodicity on each side of the interface). 
This topic has been previously studied mainly by Korotyaev via spectral theory |13l [Uj. We, on 
the other hand, use the properties of the fundamental solutions of the ID spectral problems of 
the periodic operators corresponding to each side of the interface and pose the eigenvalue problem 
as a C^-gluing problem for the decaying Floquet-Bloch solutions from either interface side. This 



approach allows us to provide some concrete conditions on Vq and the perturbation V directly 
(without conditions on the spectrum of —d'^ + Vq{x)) that ensure eigenvalue existence in the semi- 
infinite and the first finite gap of the continuous spectrum of L. Our approach is also arguably 
conceptually simpler than that of [131 E] • 

Localized waves at interfaces of two periodic (linear) structures have been also demonstrated 
experimentally in the context of electron waves in crystals by Ohno et al. |16] and for optical waves 
in photonic crystals by, e.g., Suntsov et al. jT9] . 



In detail, within the framework of the eigenvalue problem 

Li; = XiP, L = -dl + V{x), X (1.1) 
we study the following two interface problems. Firstly, an interface made of even periodic potentials 

V{x) = X{x<o}V^{x) + X{x>o}^+(^), (1-2) 
where V± has period d± > 0, i.e., V±{x + d±) = V±{x) for all x G M, and furthermore satisfies 



V±{f + x) = ViC- 
potentials 



x) for all X E [0, Secondly, an interface made of dislocated even periodic 



= X{x<o}yo{x + s) + X{x>o}Vb(a: + t), 



(1.3) 



where Vq has period d > 0, i.e., Vo{x + d) = Vo{x) for all x G M, and satisfies Vo(| + x) 



for all X € [0, ^j- The dislocation parameters are t,s € M. Here x is the characteristic function. 
Note that under the periodicity conditions the evenness of V± and Vq about x = and x = ^ 
within the periodicity cell [0,(i-t-] and [0,d], respectively, is equivalent to evenness of V± and Vq 
about X = 0. Hence, in the following we will simply require that the potentials be periodic and 
even. Unless otherwise stated, the potentials V± and Vq are continuous and hence bounded. 



One of the simplest examples of the interface ( |1.2[ ) is the additive interface 

V-{x) = Vo{x), V+{x) = Vo{x) + a, Vo{x + d) = Vo{x), Vo{-x) = Vo{x), a G M, d > 0, (1.4) 
generated by merely changing the average value of the potential on one half of the real axis. This 



example is studied in more detail in Section 3.1.1 since the conditions on eigenvalue existence 
become rather specific in this particular case. 



Schematic pictures of the two potentials (1.2) and (1.3) are displayed in Figure [1} 




Figure 1: A cartoon of example potentials V for the case (1.2) in (a) and (1.3) in (b) 



Equation (1.1) finds applications in many fields of natural science. Perhaps most notably it 



describes the wave function of an electron in a one dimensional crystal, where waves localized at 
a crystal interface are typically called Tamm states jj20j. The equation also directly applies to the 
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description of light propagating transversally to the direction of periodicity of a non-dispersive, 
lossless, linear photonic crystal which is homogeneous in the y and z directions. Suppose the 
refractive index n varies periodically in the x— direction and its mean has a jump at x = 0, such 
that n(x) = ^Jl + W{x),W{x) = -^V{x). We assume the following form of the electric field. 



^ = (0,V^(x),0)^e^('=^-'^*), 



such that the field is polarized in the y— direction, the waves propagate in the z— direction and the 
X— profile is stationary. Then Maxwell's equations exactly reduce to 



{dt - k')i^ + ^(1 + T^(x))V = 0. 



With V{x) 



-W{x) and setting A 



A; , we recover (1.1). 



Another example of an application of ( 1.1 ) is the description of matter waves in one dimensional 



Bose-Einstein condensates loaded onto an optical lattice, see Choi & Niu jJj. The density of a 
condensate is described by the wavefunction u governed by the Gross-Pitaevskii equation jJl [9l [T7] 



ihdtu 



2m 



W{x)u 



■ g\u\^u 



0, 



where, in our setting, W{x) is periodic but has a jump at x = 0. Here h is Planck's constant, m is 
the boson mass, W is the potential induced by the optical lattice and g is the scattering length. In 
the linear regime, 5 = 0, stationary waves e~'^'^^ilj{x) obey (after rescaling) equation (1.1). 

The rest of the paper is organized as follows. In Section [2] we review the needed facts on spectral 
properties of the interface-free periodic Schrodinger operators with an even potential including the 
problem of ordering of spectral band edges according to even/odd symmetry of the Bloch functions. 



Section 3.1 discusses the interface (1.2) and introduces the main tools of our analysis, namely the 
C^-matching condition and the Priifer transformation. The theory is then applied to the additive 



interface example ( 1.4 ) and numerical computations of point spectrum are performed. In Section 3.2 
we analyze the dislocation problem (1.3) for the cases s = —t and 5 = using the same tools as 
in Section [3.1 1 plus differential inequalities and variational methods. Numerical examples are, once 
again, provided. 



2 Spectrum of the Interface-Free Problem 

We review, first, some well known results on the spectrum and the eigenfunctions of the interface- 
free operator Lq := —dx + yo{x), where Vo{x + d) = Vo{x) is continuous and Vo{—x) = Vo(x). Good 
sources on the theory of the periodic Schrodinger operator are Magnus and Winkler jl5j . Eastham 
[T] and Reed & Simon tl8j. 

Lq has a purely continuous spectrum (see Theorem XIII. 90 in [18 ) consisting of bands [s2n~i, S2n] 
so that 

f^(-^o) = [J [s2n-l,S2n], 

jieN 

where s„ € M and S2n-i < S2n < S2n+i [7^- When S2n+i > S2n, we say that o"(Lo) has the 
finite gap Gn '■= (•S2n; S2n+i)- Clearly, cr(Lo) has also the semi-infinite gap Gq = (sO)Si) := 
(— oo,si). According to Floquet theory [7j the spectrum a"(Lo) can be easily found via the use of 
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Figure 2: Spectrum of Lq for Vo{x) = sin^(7rx/10). 



the monodromy matrix of the second order ODE LqV' = ^i^- Figure |2] presents the numerically 
computed spectrum of the operator Lq with Vq{x) = sin^(7rx/10). 

The ODE Loip = Xip has two linearly independent solutions, so called, Bloch functions. For real 
A 9o"(Lo) they are of the form 



(2.1) 



where A: G M \ {0} if A G int((T(Lo)) and /c G \ {0} if A G M \ ct(Lo), and ^1,2(2;; A) are real-valued 
and 2(i— periodic in x. In fact, pi^2 are either d— periodic of d— anti-periodic. If A G d{a{Lo)), the 
Bloch functions are of the form 



ijjiix; A) = pi{x; A), ip2{x; A) = p2ix; A) + xpi{x; A), 



(2.2) 



where again ^1^2(2;; A) are real and 2d— periodic in x. 

The evenness of the potential Vq{x) and the fact that only one linearly independent bounded 
Bloch function (namely ^pl{x■, X) = pi{x; X)) exists at any A G 9(c7(Lo)) imply that this solution 
must be even or odd and hence it satisfies at the boundary-points x = and x = d either Dirichlet- 
or Neumann-boundary conditions. For /c G N let {fik,Ck) denote the k-th Dirichlet eigenpair of Lq 
on [0, d] satisfying Cfc(O) = Ck{d) = and let {uf^, rjk) be the /c-th Neumann eigenpair of Lq on [0, d] 
such that ^/fc(0) = ?7^(d) = 0. The following lemma may be well known, cf. [7|, Theorem 1.3.4. 

Lemma 2.1. For the first gap edge we have si = ui. If k > 1 and if S2k S2fc+i then S2k = 
min{/_ifc, S2fc+i = max{/Xfc, 1/^+1}. Moreover, the following properties of the eigenf unctions 

are known (note that the even/ odd-property applies with respect to reflection about 





eigenvalue 


eigenfunction properties 


Dirichlet 


/^2fc-l 


even 


d— anti-periodic 


ck-i(i) = o 


Dirichlet 




odd 


d— periodic 


C2kil)=0 


Neumann 




even 


d— periodic 


^2fc-i(i) = o 


Neumann 


l^2k 


odd 


d— anti-periodic 


V2k{l) = 



Remark. Note that A G G„ can never be a Dirichlet or Neumann eigenvalue since any corresponding 
eigenfunction could be extended to a bounded solution of L^ip = Xip on M by reflection and periodic 



extension. Such nontrivial solutions cannot exist for A G G„, by (2.1) 



As we show in Sections 3.1 and 3.2, ordering between the Dirichlet and Neumann eigenvalues //^ 
and i^k+i plays an important role for existence of interface eigenvalues. It is, however, known that 
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all orderings are in general possible, i.e., for any given ordering of the Dirichlet and Neumann eigen- 
values (respecting the condition max{/xfc, z^fc+i} < min{/ifc+i, i/fc_|_2}. A; e N) a corresponding even 
potential Vq exists, see Theorem 3 in Garnett &: Trubowitz [8]. Nevertheless, the following lemma 
provides an ordering of low eigenvalues under some monotonicity assumptions on the potential Vq. 

Lemma 2.2. (a) IfV^is strictly increasing on [0, |], then V2 < Hi- The Neumann eigenfunction 
corresponding to 1/2 is strictly monotone on [0,d] and odd with respect to ^. 

(b) IfVo is strictly decreasing on [0, then fii <V2. The Dirichlet eigenfunction corresponding 
to \x\ is strictly monotone on [0, o-f^d even with respect to ^. 



The proof is based on the following result. 

Lemma 2.3. Consider a potential Vq G L°°{a, b) (not necessarily periodic, even or continuous) and 
let und be the first eigenvalue of Lq = —d"^ + Vo on [a,b] with the boundary condition u'{a) = = 
u{b), whereas kdn denotes the first eigenvalue of the same differential operator but with boundary 
conditions u{a) = = u'{b). Then 

min{KjvD) i^dn} = min i^J v''^ + Vo{x)v'^ dx : v £ H^{a, b) has a zero and J dx = 1^ . (2.3) 

Moreover, ifVo is strictly increasing on [a,b] then K]\j£) < hdn o-nd any eigenfunction for with 
u{a) > is strictly decreasing on [a,b]. If Vq is strictly decreasing on [a,b] then kdn < i^nd cLnd 
any eigenfunction for kdn with u'{a) > is strictly increasing on [a,b]. 



Proof. The proof is inspired by a similar result in Bandle et al. [2j. Note first that the set, on 
which the minimization is performed, is weakly closed in H^{a,b) due to the compact embedding 
H^{a,b) C[a, 6]. Hence a minimizer of the right-hand side of (2.3) exists. We denote it by U. 
Let us also denote the value of the minimum by k. The proof is now divided into five steps: 

Step 1: U has exactly one zero on [a,b]. Since U possesses at least one zero G we have 

U G = {v £ H^{a, b) : v{xo) = 0}. Clearly U is then the minimizer of 



min < / v'"^ + Vq{x)v'^ dx : V £ HxQ, I v"^ dx = 

(.J a J a 

and therefore U satisfies the Euler-Lagrange equation 

- U" + Vo{x)U = kU in (a, xq) U (xq, b) (2.4) 

with boundary condition 

U'{a) = U{xo) = U'ib) = (2.5) 
where in case xq G {a, 6} one of the two Neumann conditions is dropped. Note that 

j-b pb 

/ U'v' + Voix)Uvdx = K Uvdx for all v G i^^o- (2.6) 

J a J a 

Now assume for contradiction that U has a second zero xi 7^ xq. Then (2.6) holds also for all 
V G and since H^{a,b) = H^^ © Hx2i we find that (2.6) holds for all v G H^{a,b), i.e., U is a 
Neumann-eigenfunction. The same applies for \U\, which is also a minimizer of (|2.3[). But then U 
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must be the first Neumann-eigenfunction of Lq on (a, b) and it therefore has no zero on [a, b]. This 
contradiction shows that U has exactly one zero in [a, b]. 

Step 2: K is strictly less than the second Neumann- eigenvalue V2 on [a, b] . Since the second Neumann 
eigenfunction 772 has one zero in [a, 6], we find k <V2. Suppose for contradiction that n = V2- Testing 
the equation for 772 with 7/^ = max{r/2,0} we obtain 



(r?2+ Y + Vo{x){rit) dx = U2 {4) dx 



and thus ^ minimizer for (2.3) and must have a unique zero by Step 1. However, clearly 7/2 

has a continuum of zeros. Therefore we can conclude that k < V2. 

Step 3: U has its unique zero either at x = a or at x = b. If we suppose for contradiction that the 



unique zero xq lies in the open interval {a,b), then we obtain the Euler-Lagrange equation (2.4) 



with boundary condition (2.5 1 . By rescaling the minimizer U suitably on [a, xq] we can achieve that 



the rescaled function C/ is a C^-function on [a, b] solving the equation (2.4) pointwise a.e. on (a, b). 
Hence, the rescaled function [/ is a Neumann-eigenfunction with one interior zero, i.e., k = 1^2 in 
contradiction to Step 2. 

Now the claim of the lemma about the value of the minimum is immediate. 

Step 4- ordering of k^d, kdn- We are using the following rearrangement result of Hardy, Little- 



wood, Polya [To]. Let v,w be non-negative and measurable on [a,b]. If 



are the increasing 



rearrangements of v,w, then J^vwdx < v'^w" dx. Moreover, if v is strictly increasing, then 
equality holds if and only if w = wK A similar statement holds for the decreasing rearrangements 
v*,w*. Note, that the non-negativity of v,w can be replaced by boundedness. 

A simple corollary of the Hardy, Littlewood, Polya inequality is the following: suppose 
is strictly increasing and both V and w are bounded. Then 



/ Vw*dx< I 

J a J a 



Vw dx 



(2.7) 



with equality if and only if w = w* . The proof follows immediately from the observation that 
Let Vq be strictly increasing on [a,b]. Suppose for contradiction that kjjn < i^nd and let U 



be an eigenfunction corresponding to kd^, which by (2.3) is also a minimizer of the variational 



problem in (2.3). We may assume U to be non-negative, since \U\ is also a minimizer of the 



corresponding variational problem and K£)]\i is a simple eigenvalue. Let now U* be the decreasing 
rearrangement of U on [a,b] and note that (C/^)* = (U*)'^. Since for the decreasing rearrangement 
we have fl^{U*'Ydx < fl^{U'fdx, cf. Kawohl HI], we obtain by apphed to Vq and C/^ the 

relations 



{U*ydx 



dx = 1, 



{u*'y + Vo{x){u*ydx< / {u'y + Vo{x)u'dx. (2 



Therefore U* , which satisfies U*{b) = 0, is also a minimizer of 



in (2.8). But since Vq is strictly increasing, the sharp form of (2.7 



2.3) and hence equality has to hold 
implies that U = U* which by 



U{a) = implies the contradiction that U must be identically zero. Hence katd < i^dn- Moreover, 
(|2.8[) shows that any non-negative minimizer U for ktvd satisfies U = U*, i.e., U is decreasing, and 
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by using the differential equation for U and tlie strict monotonicity of Vq it is easy to see that in 
fact U is strictly decreasing. 

If Vo is strictly decreasing on [a, b] then a similar argument based on replacing U by its increasing 
rearrangement shows that hdn < f^ND- D 



Proof of Lemma 2. 2 



[0, is the eigenfunction for kdn of Lemma 



eigenfunction for knd 
from Lemma 12.31 



Consider the Dirichlet-e igen function 
)n for K£)j\f of Lemma 2.3 
Hence fj-i = k^n and V2 



By Lemma 2.1 its restriction to 



Likewise, the restriction of r/2 to [0, is the 
und- The statements (a) and (b) then follow 

□ 



3 Interface Problems 



Let L be the operator in ( 1.1 ) defined on the dense subset (M) of L (M). We investigate next the 



existence of eigenvalues of L for the interface potentials ( 1.2 ) and ( |1.3[ ). These examples fall into a 
larger class of potentials, namely V{x) = X{x<o}^i{x) + X{x>o}^2{x) ^ where Vi^2{x + di^2) = ^1,2(2;) 
for some di,2 > but where Vi^2 may not be even in x. Clearly, all solutions of {—d'^ + V{x))'ijj = \ip 
are then 

where ■i/'i are Bloch functions of (—9^ + Vi^2{x))il) = A^/^, respectively. As decaying Bloch functions 
V'i exist only in spectral gaps of —9^ + 14,2(2;), respectively, eigenvalues of L can exist only within 
intersections of the gaps of o"(— 9^ + Vi{x)) and a{—d'^ + V2{x)). Note the following additional 
information on the spectrum of L, which for our purpose plays no further role: the essential 
spectrum of L is the union of the essential spectra of —d"^ + V\{x) and —d^ + V2(x), cf. Korotyaev 
[T?j . As a result, no embedded eigenvalues of L exist. 



3.1 Point Spectrum for Interfaces Made of Even Potentials 



The eigenvalue problem (1.1) with (1.2) can be viewed as the system 



-dl^i + V-{x)ii = Xip for X < 0, 
+ V'-^-{x)^|J = Xijj for x > 



coupled by the C^-matching conditions 

iP^O-) = ^(0+) and V'(O-) = V''(0+). 
As stated in Section [T] the functions V±{x) are continuous, even and d-t -periodic. 



(3.1) 



(3.2) 



Based on the knowledge of the fundamental solutions in (2.1), (2.2) we conclude that an L - 



integrable solution of (1.1) with (1.2) can only exist if A lies in the intersection of the resolvent 



sets, i.e., in the intersection of the spectral gaps of L_ and i.e., if A G H for some 
n, m G N U {0}, where is the n-th spectral gap of L± respectively. 

For A € G^ n G~ with some n,m >0 any localized eigenfunction ip of L, therefore, has to be of 
the form 

'iPix; A) = x{x<o}i^~{x; A) + x{x>o}V'+(a;; A), 
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where 

V'±(^;A)=p±(x;A)eT«W^ (3.3) 
with k{\) > and p±{x; A) being 2 d-i-— periodic in x. The functions p± are restrictions of either pi 



or p2 in (2.1) with Vq = V± to the half-Une M-t respectively. 



An important remark is that, due to the linearity of the problem, the matching conditions (3.2) 
together with an appropriate scaling are equivalent to 

fl+(A)=i?_(A), where R^{X) = "^j^^ (3.4) 

V'±(0; A) 

and the prime denotes differentiation in x. 



We determine existence of solutions to (3.4) via the intermediate value theorem and by mono- 



tonicity of the functions i?±(A). The monotonicity then also implies uniqueness. 

Lemma 3.1. Within each gap and G^,n > 0, the functions Rj^ and are continuous 
functions of X , which are strictly increasing and decreasing respectively. 

Proof. Let us start with the proof for i?_(_(A). Under the Priifer transformation, cf. Coddington &; 
Levinson [5] 

ip+{x; A) = p{x; A) sm{6{x; A)), tp+{x; A) = p{x; A) cos(^(x; A)), 

the equation L^ijj^ = Xip^ becomes 

e' = l + {X-V+{x) -l)sm'^{e), 

p' = -p{X - V+{x) - 1) sin(^) cos(0), 

where the prime denotes differentiation in x. Clearly, 6 and p are continuous functions of both 
variables a; G M and A G and since i?+(A) = cot(0(O; A)), the function -R+(A) is continuous in A 
provided V'+(0; A) has no zero in the interior of G+. Note that if ■0+(O; A) = 0, then by evenness 
of V+ and the reflection symmetry of the problem = A'i/'+, the solution iIj+{x] A) defined in 



(3.3) on X > could be extended to a solution on x G M via ^/'+(— x; A) = — V'+(x; A). This solution 
would decay exponentially at both infinities and A would, thus, be an eigenvalue of L-(_, which is 
impossible. Hence continuity of R+{X) is proven. 



Now let us prove the monotonicity. Due to the form of '0+ j see ( |3.3| ) , we have 

p(2d+) = ^(V.+ (2d+))2 + (V.;(2d+))2 = e-2'^+XO). (3.5) 

Define now z(a;) := ||(x;A). The function z satisfies z' = z(A — V+(x) — 1)2 sin(^) cos(0) +sin^(^) = 
-2^z + sin2(6i). Therefore, 



Because cot{9) = and due to the periodicity ^^(^^2d+'A) ~ ^+^(x'-X] ^^ have 6{2d^; A) = ^(0; A) + 
rmr, where due to continuity the value m G Z is independent of AP] Hence, z{2d^; X) = z(0;A). 



Using (3.5) and (3.6), we thus obtain 



^In fact, it can be easily seen from Sturm oscillation theorem that m = 2n, where n is the index of the gap G^. 
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Because k > 0, we get z(0) < and conclude that 6{0; A) is strictly decreasing throughout 
Therefore, R+{\) = cot(0(O; A)) is strictly increasing with respect to A throughout G+. 



In order to prove strict monotonicity of R-{X), note that (3.5) is replaced by p{—2d-) = 
g-2d_K^^Q^ and in (3.7) the value 2d+ is replaced by —2d- both in the arguments of the func- 



tions z and p and in the upper limit of the integral. This leads to the conclusion z(0) > which 
means that i?_(A) is strictly decreasing with respect to A. □ 

In order to apply the intermediate value theorem and prove crossing of the graphs of R+{X) and 
i?_(A), we use their continuity within each gap and their limits as A approaches a gap edge. 

Lemma 3.2. Let s G {si,,S2,...} be one of the boundary-points of the spectral gaps of L± 
respectively. If s corresponds to a Dirichlet-eigenvalue of L± on [0,d], then lim^^^^^^Q± 
|i?_l_(s)| = oo respectively, and if s corresponds to a Neumann-eigenvalue of L± on [0,d\ then 
lim^^^ agG* ^±(^) = R±{s) = respectively. 



Proof. We only consider the "+" case. Let A^ G A^ ^ s be a given sequence. Due to (3.3) the 
functions '(/'+(•, Afc) have the form 



ip+ix; Afc) = p+(x; Afc)e" 



where w.l.o.g. we may assume Afe)||Loo = 1, which implies ^^fc) IIl°°([o,oo)) ^ 1- On 

every compact subinterval [0, b] C [0, cxd) the H^-norm of tp+{-, Xk) is uniformly bounded in k and 
hence along a subsequence (again denoted by Afc) the functions ip+{-; \k) converge in H^{[0,b]) 



(and hence, by the differential equation (3.1) also in H {[0,b])) to a solution v of L+w = sv with 
II^I|l°°([o,6]) — 1- Since this holds for every 6 > 0, the function f is a bounded solution of L-^v = sv 



on [0,oo) and therefore coincides with the bounded periodic Bloch function pi{x;s) in (2.2). The 



convergence of i?+(Afc) is now obvious by the embedding H'^{[0, b]) into (^""^([O, b]). □ 

To make the picture of the behavior of R± complete, it remains to determine their behavior at 
the lower end of the semi-infinite gap Gq , i.e. as A ^ — oo. 

Lemma 3.3. Let V± be bounded potentials (not necessarily even, periodic or continuous). Then 
R±{X) =poo as X ^ —oo. 



Proof. The proof is, as for Lemma 3.1 shown only for i?+ with the one for i?_ being completely 
analogous. We rescale the Bloch function tp-^{x; A) so that tp+{0; A) = 1. Note that this is possible 
if and only if il^+(0; A) ^ 0, which we show to be true for all A < inf V+. Suppose that il^+(0; A) = 0. 
Testing (L_|_ — Xjijj-f- = with ip+ over x £ [0, oo), we get 

/ {iJ+fdx+ / {V+ - X)ipldx = 
Jo Jo 

and, therefore, A > inf V+. 

Let now A = —u^ for some u > 0, s.t. — i^^ G G^ and — z/^ < inf V+, and define 

(/>^(a;) := V+(x;-i^^) -e-^^ (3.8) 

We have 

£ = u^(pu + V+i;+, 0.(0) = 0. (3.9) 
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Since -R+(A) = z^^) = '0+(O; —v'^) = — z^ + i;^'j^(0), we need to determine the behavior of (l)'^{0) 

as v ^ oo. Using the Green's function, we solve (|3.9l) to obtain 



sinh(z^t)F+(t)V'+(t;-z^ )dt + sinh(z^a;) / e~''^V+{t)ij+{t; -iy^)dt 



Therefore, (p'^{0) = - e-^*F+(t)V'+(t; -i^^) dt and 

\Km < ||y+||L-||e-'^-||i2(o,oo)IIV'+(-;-i^' 



\v+ 



L2(o,oo) 



Il2(0,oo)- 



(3.10) 



In order to estimate HV'+Cs )IIl2(o,oo)) (3.9) yields 



^^II'^^IIl2(0,oo) 



^'1^111,2(0,00) 



V+{x)ip+{x; -v'^)(t)u{x)dx, 



implying 1101/11^2 (0,00) ^ ll^+lli°°llV'+(-; -i^^)I|l2(o,oo)I|</'i/||l2(o,oo) 



and 



^i^llL2(0,oo) < ^II^+I|l°°II'0+(-; -i^^)llL2(o,oo)- 



Therefore ^ and together give ||V'+(-; -y^)\\L^(^o,oo) < ^ + ^11^+11 

If > 1 1 V+ 1 1 2, 00, we have the estimate 



^+(-;-^-^)iil2(o,oo)< II 



(3.11) 



Il2(0,oo)- 



(3.12) 



Finally, combining (3.12) and (3.10), we arrive at the bound 

i2u)-^\V. 



L°° 



l-u-^\\V+ 



\L°° 



which implies i?+(— i^^) = —v + (j)'^{0) — > —00 as — > 00. 



□ 



The behavior of the ratio functions R±{X) for the two examples = V- = sin^(7r3;/10) and 
V+ = y_ = cos^(7rx/10) is summarized in Figure [sj Note that Lemmas 2.2 3.1 3.2 and 3.3 imply 
the behavior only for A < S3. The rest in Figure [3] is obtained without a rigorous proof from 
numerical computations of the gap edge eigenfunctions. 

By the intermediate value theorem and based on the behavior of R±, we now obtain the following 
theorem, which has already been observed by Korotyaev |14] . 

Theorem 3.4. Let G^,G^ be two gaps in the spectrum of L„ and respectively, such that 
n 7^ 0. Then the following two statements are equivalent: 

(a) 3 A G G^ n such that A is an eigenvalue of L. 

(b) Either G~ = {fJ-n , ''n+i) ^ G+ = (z^++i,^+) or G" = (z/-_^i, //"), G+ = (/x+,f++i). 



In the affirmative case the eigenvalue is also unique. 



10 



(a) 







1 1 

G 




1 

G, 
1 


1 1 

G 




=3 

h% 






-00 







CO -co 




/^oo 




Rp.) 


00 


\l 
1 1 


1 


-00 CO 
1 1 




0\-oo 

1 







0.2 0.4 


0.6 


0.8 , 1 
A 


1,2 


1.4 


1 










(b) 









,6 



T 



-CO -CO -co -oo/f 

o3\^0co \^ Oco\^0 CO \^ 

I I I I I I I 

0.2 0.4 0.6 0.8 1 1,2 1.4 1.6 

7. 



Figure 3: Behavior of the ratio functions R±{\) within gaps of a{L±) for V+ = V- = sin2(^j;/10) 
in (a) and for V+ = y_ = cos^(7rx/10) in (b). The arrows denote the monotonicity type. 





Figure 4: A cartoon of the graphs of ii_(A) and i?+(A) when the former one of the two conditions 
in Theorem 3.4 holds, (a) No solution to (3.4) without monotonicity of i?_ (hypothetical case), 
(b) Existence and uniqueness of the solution to (3.4) on any H with monotonicity of R±. 



Remark. Note that besides continuity and the limit values of R±{X) their monotonicity is also 
needed to fulfill the conditions of the intermediate value theorem. Without monotonicity the ranges 
of the functions R+{X) and i?_(A) on the intersection G~ n could be completely distinct, see 
Figure|4](a). With monotonicity of i?-|-(A) we, of course, obtain also uniqueness of solutions to (3.4). 

Let us call the gap (fJ^t^'^n+i) ^ DN-gap and the gap {i^^+i^ l^n) ND-gap. The semi-infinite 



gap belongs to the class of DN-gaps. The existence part of Theorem |3.4| can then be formulated as 
follows: 

Whenever a DN/ND-gap of intersects an ND/DN-gap of L^, respectively, a unique 
eigenvalue of L exists in this intersection. 
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3.1.1 Example: additive interface 



The additive interface problem (1.1) with (1.4) is equivalent to (3.1) with L_ 
and La 



Lo ■.= -dl + VQ{x) 



Lq + a. 

Because a{LQ + a) = (t(Lo) + a, we have = G~ + a, and because the Bloch functions of Lq 
at the spectral parameter A are the same as the Bloch functions oi Lq + a at A + a, to check the 



conditions of Theorem 3.4 one only needs to know cr(Lo) and symmetries (even/odd) of the Bloch 



functions of Lq at the gap edges s„- 



The existence part of Theorem 3.4 can now be formulated as follows: 



Whenever a shifts the spectrum of Lq so that a shifted DN/ND-gap intersects an (un- 
shifted) ND/DN-gap, respectively, a unique eigenvalue of L exists in this intersection. 



Theorem 3.4 has several interesting and rather specific corollaries for the additive interface case. 



Firstly, clearly, if |a| < q*, where a* := inf„gN(s2n — S2n-i) stands for the width of the narrowest 
spectral band of Lq, the shift a is too small to make even the two gaps lying closest to each other 
overlap. 



Corollary 3.5. // |a| < a* := inf„gN(s2n — 



:n-lj 



then L has no eigenvalues. 



In the rest of this section Gn denotes the n-th spectral gap of Lq. As Lemma 2.2 dictates, when 
Vq is strictly increasing on [0,(i/2], the first finite gap Gi = (s2,'S3) is an ND-gap and thus if a 
shifts the semi-infinite (DN) gap Go so that Gq + a intersects Gi, an eigenvalue exists. Obviously, 
the infimal value of a > achieving such an intersection is the width of the first spectral band 
■S2 — si. Since Go is semi-infinite, there is no upper bound on a and if a > S2 — si, the intersection 
is always nonempty. On the other hand, when Vq is decreasing on [0,(i/2], Gi is a DN-gap and 
the intersection of Go + a and Gi contains no eigenvalues. As the next Corollary clarifies, for 
a < — (s2 — si) the situation is similar. 

Corollary 3.6. Let Vq be strictly increasing/strictly decreasing on [0, |]. // \a\ > S2 — si then a 
unique eigenvalue/no eigenvalue of L exists in Gi n (Go -|- a) for q > and in Gq n (Gi -|- a) for 
a <0. 



Remark. For the case of the additive interface it is possible to show that the number of eigenvalues 
of L is finite for any a G M based on the asymptotic behavior of gap locations and gap widths. 
Indeed, based on Theorem 4.2.2 in [7] the center of the n-th gap behaves like an? + o{n) as n ^ oo 
with the constant c G M dependent on Vq. The gap widths, on the other hand, tend to since 
they build an sequence, see Theorem 3 in [8]. Therefore, asymptotically, the n— th gap has the 
form cn? + J„, where both inf(J„) and sup(J„) behave like o{n). For a given a G M infinitely many 
eigenvalues are thus possible only if for infinitely many pairs (m, n) G N x N with n ^ m there exist 
Sn S Jn and tm G Jm such that 



cn -|- Sr. 



cm + tm + Oi. 



(3.13) 



As for 



n 



m no eigenvalues exist, we can rewrite (3.13) as 

/ , \ Q. -\- tjYi Sn 

c[n + m) = 



n — m 

Clearly, the right hand side is o(n -|- m) while the left hand side is not. Thus only finitely many 



solutions of (3.13) exist. 
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For general interface problems (with V{x) = X{x<o}^-{^) + X{x>o}^+(^)) the question of finite- 
ness of the number of eigenvalues seems open. Due to Theorem 3 in |i8j there are, for example, po- 
tentials V- and V+ with equal gap lengths and opposite DN/ND 'polarities'. //, in addition, the lo- 
cations of the gap centers were identical, there would be an eigenvalue in each gap G~ = G^, n G N. 
However, it seems to be an open problem whether such potentials V- and V+ exist. 



Numerical results The point spectrum of the additive interface problem with the potential 
Vo{x) = sin^(7ra;/10) has been computed using a 4th order centered finite difference discretization. 
The eigenvalues are plotted in Figure [5] for a range of values of a. The shaded regions are the union 



of spectral bands of Lq and Lq -|- a. The results agree with Theorem 3.4 




Figure 5: Numerically computed point spectrum of L with Vq{x) = sin^(7rx/10) for a range of 
values of a. The union of spectral bands of Lq and Lq -|- a is shaded. The inset blows up the region 
near A = si, a = 0. Eigenfunctions for the labeled points are plotted in Figure |6j 



In Figure |6] we plot eigenfunctions corresponding to nine selected eigenvalues in Figure |5] Note 
that the decay rate of the eigenfunctions is often very different on either side of the origin. 

For the potential Vo(x) = cos^(7rx/10) it is clear from the numerically obtained Figure[3](b) that 
the intersections Gj R {Gk + a), j, k £ {0, ... ,3} contain no eigenvalues because the gaps Go and 
Gi are DN-gaps and so seem to be G2 and G3. In other words, based on the numerics, the additive 
interface problem (1.1), (1.4) with Vo(x) = cos^(7rx/10) has no eigenvalues on (—00,83]. Note that 
our analysis guarantees non-existence of eigenvalues in ( — 00,54]. 



3.2 Point Spectrum for Interface Problems Made of Dislocated Even Potentials 



For the dislocation interface (1.3) we restrict our attention to the two representative cases t = —s 
and s = 0. 
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Figure 6: Eigenfunctions corresponding to the 9 labeled eigenvalues in Figure [5] 



3.2.1 Symmetric Dislocations 



Here we study the eigenvalue problem (1.1) with (1.3) in the case where t = —s, t € (0,d). This 
can be done via the system 



:= -<9^V* + Voix - t)ip^ = \ip^ for x < 0, 
L\^t ■= -dlijj* + Vq{x + t)^* = A^* for X > 



coupled by the the C^-matching conditions 

^*(0-) = V*(0+) and 



dx 



dx 



V'*(o+). 



(3.14) 



(3.15) 



First note that the spectrum a{Lt) of the operator Lt := —9^ + Vo(x + i) on M is identical to the 
spectrum (j{Lq) of Lq = —d"^ + Vq on M and we have = G~. Moreover, the Bloch functions 
Tpi 2 of Lt for A G M \ (y{Lt) = M \ (t{Lq) are just shifts of the Bloch functions of Lq, i.e., ip\{x] A) = 



■0^(x + t; A), i = 1,2. Therefore, an L^-solution of ( 1.1 ) with (1.3 ) can only exist if A ^ cj(Lo). For 



such A any localized eigenfunction ip^ of (1.1) with (1.3) must take the form 



'4^\x] A) = X{x<o}'4'-{x] A) + X{x>o}V'+(^; -^)' 
where ^^(x; A) are those Bloch functions of L±t^ which decay on M^, respectively. 
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As in Section [3.11 we introduce the ratio functions 



R±{x; A) 



so that the matching conditions (3.15) are equivalent to i?'^(0; A) = R^_{^; A). Due to the fact, that 
the Bloch functions are just shifts of the Bloch functions we see that i?+(x; A) = i?1_(a;+t; A) 
and R^_{x; A) = W_{x — t]\). Thus, the matching condition (3.15) amounts to 

R%(t-\) = ii° (-t;A). 

Finally, the evenness of the potential Vq and the fact that only one linearly independent Bloch 
function decaying at +oo exists, imply that 'iIj^_^{x] \) = ±^^1(— x;A) since A ^ cr{LQ), and hence 
R^{t; A) = —R^(—t; A) so that finding an eigenvalue of (1.1) with (1.3) amounts to finding a zero 
or a pole of R^{t; X) for some t G {0,d). This is done below via the intermediate value theorem 
and monotonicity properties of the function R\{t; A). 

For simplicity we write in the following R{t; A) instead of i2+(t; A). First, we need to generalize 
Lemma 3.1 on the monotonicity and continuity of R{t; A) or the corresponding Priifer angle 6{t; A) 
as a function of t and A. Suppose ip £ -L^(0,oo) solves L^tp = Xip. We apply again the Priifer 
transformation given by 

ip{x; A) = p{x; A) sin(^(x; A)), ip'{x; A) = p{x; A) cos(^(x; A)), 

which transforms the equation Lqi^ = Xip into the system 

e' = l + {X-Vo{x) -l)sm'^{9), (3.16) 
p' = -p{X - Vq{x) - 1) sin(0) cos(^), (3.17) 



where the prime denotes differentiation in x. Note that ( |2.1[ ) implies 2d-periodicity in t of i?(t; A). 
Hence 9{t + 2d; A) = 6{t; A) +rmr, where m is an integer which is constant in A within each spectral 
gap. In fact, it can be shown by the Sturm oscillation theorem that m = 2n when A G G„. 



In the subsequent arguments we use the following result on differential inequalities, cf. Wal- 
ter [22], which we quote in a slightly simplified way. Functions v,w satisfying (3.18) below are 
called sub-, supersolutions, respectively. 

Lemma 3.7. Let f : [a,b] x R ^ M be continuous and continuously differentiahle with respect to 
the second variable. Ifv,w£ C^[a,b] satisfy 

v' < fitjv), w' > fitjw) on [a,b] with v{a) < w{a), (3.18) 

then V < w in [a, b]. More precisely, either v < w in (a, b] or there exists c G (a, b] such that v = w 
on [a, c] and v < w on (c, b] . Moreover, if one of the differential inequalities holds strictly almost 
everywhere in [a, b], then v{t) < w{t) holds for all t £ {a, b]. 

Lemma 3.8 (Monotonicity in A). LetGn = {s2n, S2n+i),n' >0 be a fixed gap. For (t. A) G [0,d]xG„ 
the function R is continuous except in the set S = Ute[od] '^t' '"^^^'^s for each t either St = 9 or 
St = {{t, At)} and 

lim R{t; X) = +oo, lim R{t; X) = -oo. (3.19) 

A— ►At— A^At + 

For a fixed t the function R(t; A) is strictly increasing for X G Gn if St = 0, and strictly increasing 
for X G {s2n, Xt) and for X G (A^, S2n+i) if St = {{t, Xt)}. Moreover, if X, p £ Gn and St / 0, then 
X < Xt < p implies R{t; X) > R(t;p). Consequently, for all t G [0,d] we have that X ^ p implies 
R{t;X)^R{t;p). 
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Proof. As we have seen in Lemma |3.1[ the Priifer- variables 9 and p are continuous functions of 
both i E M and A E Gn- Since R{t; A) = cot{9{t; A)), the function R(t; A) is continuous except for 
those values, where 6{t; A) passes through kir, k £ Z. Since R{t; A) is strictly increasing in A at 



points of continuity by Lemma 3.1, the relation (3.19) follows. The fact that there is at most one 



blow-up point Xtg with respect to A will follow from the next statement. Let Aj^ be a pole and 
X < XtQ < fJ- and suppose for contradiction that R{tQ; A) < R{to; p). By lowering p if necessary and 
keeping the order A < Atg < /U, we may achieve R{tQ;X) = R{tQ;fi), i.e., there exists k £ Z such 
that 9{to; A) = e{to; fi) + kn. Note that 

e'it; A) = 1 + (A - Voit) - 1) sm\0{t; A)), 

0'{t; ^) = 1 + - Vo{t) - 1) sm\0{t; ^)) > 1 + (A - Voit) - 1) sin2(^(t; p)) 



for almost all t > Iq. By the comparison principle of Lemma 3.7 we obtain 9{t; A) < 9{t; /i) + k-ir 
for all t > tQ. Here we have used that 9 and 9 + kn solve the same differential equation. It follows 
in particular, that 

9{to; A) + mvr = 9{to + 2d; A) < 9{to + 2d; fi) + kfr = 9{to; fi) + {m + k)TT 

contradictory to our assumption 9{to; A) = 9{to; n) + kir. This proves the lemma. □ 

Corollary 3.9. For t = —s the number of dislocation eigenvalues in any gap Gn,n > 0, is 0, 1 or 
2. If there are 2 eigenvalues, then one of them has an even and the other one an odd eigenfunction. 



Proof. It follows from Lemma 3.8 that for fixed t the function R(t; A) as a function of A can have 
at most one zero and at most one pole. □ 

Lemma 3.10 (Monotonicity in t). Suppose Vq is an even, d — periodic -function. Let G^ — 
is2n, S2n+i),n' > 0, 6e a fixed gap and let X E dGn- 



(a) If Vq is strictly increasing on [0, then either 9{t; A) is strictly increasing for t E [0, d\ or 
there exists to E (0, ^) such that 9{t;X) is strictly increasing for t E [0, to] U [d — to,d] and 
strictly decreasing for t £ [to, d — to]. 

(h) If Vo is strictly decreasing on [0, j]; then either 9{t; A) is strictly increasing for t E [0, d\ or 
there exists to E (0, such that 9{t;X) is strictly decreasing for t E [0,to] U [d — to,d] and 
strictly increasing for t £ [to, d — to]. 



Note that in both cases, 9{t;X) can change monotonicity with respect to t only once on [0,d/2]. 



Proof. We gi ve th e proof in case (a). The proof for case (b) needs only minor modifications. Recall 
that for A E dGn the evenness of Vq implies 9{^; A) = /c| for some k £ Z. Hence 



from Lemma 
we have 



2.1 



d 



+ s;X 



kTr-9 



s; X] Vs E 



0, 



d' 



since both sides satisfy the differential equation (3.16) with Vo{x) = Vo{d/2 + s) = Vo{d/2 — s), and 



have the same initial values at s = 0. In particular 

'd 



s; X\ Vs E 



0, 



(3.20) 
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due to the evenness of Vb(x) about x = d/2 (implied by d— periodicity and evenness about x = 0). 
In any of the two cases, the monotonicity of in [0,i] has its counterpart in [|, c?]. Differentiation 



of (3.16) with respect to t yields 

{ey = 2(A 



Vo{t) - 1) sm{e) cos{e)e' - V^it) sia^ie). 



If e'{to) > for some to G (0,d/2], then by ( |3.20[ ) also e'{d - to) > 0, so that Lemma [3^ applied 
to t> := and w := 6' on [d — to, d] (note that V^' < a.e. on [d/2, d]) implies 6*' > on (d — to, d] 



and by (3.20) also on [0,to). Below we show that such to exists. Let to be chosen maximal with 
these properties. If to = d/2, then 6 is strictly increasing on [0,d]. If to < d/2, then 0'{to) = 
and Lemma 3.7 applied to v := 6' and u; := on [to, d/2] (note that V^' > a.e. on [0,(i/2]) gives 
61' < on (to, d/2] and by ( |3.20D also on [d/2,d - to). Consequently, 6' > on [0, to) U{d-to,d], 
and 6' < on (to, d — to)- 

Note that the case 6' < throughout (0, d) is impossible since then 6{0) ^ 9{d) and thus 6{s) 
is a multiple of n for some s S [0, d] (because ^(0) and 9{d) are multiples of | by the remarks 



before Lemma 2.1), whence 0' > in some neighborhood of s due to the differential equation for 
6, contradicting 9' < on (0,d). This proves the lemma. □ 



Suppose Vo satisfies the basic assumptions, i.e., it is even, d-periodic and contin- 

-oo,si). 



Theorem 3.11 

uous. Let s = —t in (1.3) and consider the semi-infinite gap Go 



(a) If Vq is strictly increasing on [0,d/2], then there is no/exactly one dislocation eigenvalue in 
Go fort £ [0,d/2] / {d/2,d) respectively. 

(b) If Vo is strictly decreasing on [0,d/2], then there is exactly one/no dislocation eigenvalue in 
Go for t G (0,d/2) / [d/2,d] respectively. 

Proof. It suffices to prove part (a), since (b) follows from (a) via shifting the potential by the 
half-period | due to the evenness of Vo{x) about x = d/2. Recall that the first band edge si 
is a Neumann eigenvalue. The first Neumann eigenfunction u is positive, and hence, due to the 
d— periodicity and evenness of Vo it has an extremum at x = d/2. It can be thus viewed as the first 
Neumann eigenfunction on the interval x G [0,d/2], i.e., the minimizer of the energy 



d/2 



v'^ + Vo{x)v^ dx. 



rd/2 

where w e //^(O, d/2) with / v^dx = l. 

Jo 



As the decreasing rearrangement u* of u decreases the energy, u has to be decreasing, i.e., u'{x) < 0, 
on [0,d/2]. In fact, u'{x) < on (0,d/2). If n'(^) = for some ^ G (0,d/2), then due to positivity 
of u the function R satisfies i?(0; si) = i?(^; si) = R{d/2; si) = 0, hence R, and in turn 9, change 
monotonicity at least three times on (0,d/2), which is impossible by Lemma 3.10 Therefore 
R{t;si) < for t G {0,d/2), R{t;si) > for t G {d/2,d), and R{t;si) 

Ri{0;X) - 



Recall now from Lemma 



3.3 



for t G {0,d/2,d}. 
cxD as A ^ — oo for any t G [0,d]. 



that R{t;X) 

Moreover, R{t; A) is continuous in A G Go because continuity can be broken only by a pole. But 
because R{t; A) — > — oo as A — > — oo and i?(t; A) is increasing in A within each continuity segment, 
a pole would mean that i?(t; A) takes the same value for some Ai 7^ A2 G Go, which is impossible 
by Lemma |3.8| 



As a result R{t; X) stays negative for t G {0,d/2) throughout A G Go, goes through once for 
t G {d/2, d), and takes the zero value at A = si ^ Go for t = d/2. □ 
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Theorem 3.12. Suppose Vq is an even, d-periodic -function, let s = —t in (1.3) and consider 
the first finite gap Gi = (s2, S3). 

(a) Suppose Vq is strictly increasing on [0,1]- Then Gi = {i'2^l^i), CLnd the second Neumann- 
eigenfunction is strictly monotone on [0,d]. For the first Dirichlet-eigenfunction u we have 
the alternative: 

(al) u is strictly monotone on [0, Then there is exactly one dislocation- eigenvalue in Gi 
for t G (0, d)\{^} and none for t = ^. 

(a2) u changes monotonicity on [0, |] exactly once at the extremal point do G (0, |). Then 
the number of dislocation- eigenvalues in Gi is as follows: 



dislocation parameter 


t G (0,4) 


te[do,^], 


t G {i,d-do) 


t e[d- do,d) 


t = d 


number of eigenvalues 


1 





2 


1 






(b) Suppose Vq is strictly decreasing on [0, j]. Then Gi = (/ii,f2), and the first Dirichlet-eigen- 
function is strictly monotone on [0, 1]- For the second Neumann- eigenfunction u we have the 
alternative: 

(bl) u is strictly monotone on [0, g]- Then there is exactly one dislocation-eigenvalue in Gi 
for t G (0, d)\{^} and none fort = ^. 

(b2) u changes monotonicity on [0, exactly once at the extremal point do G (0, |). Then 
the number of dislocation- eigenvalues in Gi is as follows: 



dislocation parameter 


t G (0,4) 


te[do,1), 


t-^ 
^ 2 


t G (f,c?-4) 


t G [d-dQ,d\ 


number of eigenvalues 


2 


1 





1 






Proof. As in Theorem 3.11| it suffices to prove part (a) when, in addition, the roles of V2 and m are 



switched in the proof of (b). The strict monotonicity of the second Neumann eigenfunction and 



the fact that Gi = (z^2)/^i) was already stated in Lemma 2.2 For the monotonicity alternative of 



the first Dirichlet eigenfunction u (which can be assumed positive on (0,(i)) recall that 

u{x) = p{x; ^i) sin 9{x; /ii), u\x) = p(x; /ii) cos 9{x; fii). 

We can assume that 9{0; fxi) = 0, 0(|;^i) = | and that 6{x;fii) ranges in [0,7r) for x G [0, |]. 
According to the monotonicity alternative for 6 in Lemma |3.10| there are two possibilities: either 
6{x; jjLi) is increasing and hence cos{9{x; /ii)) > for x G (0, g), or 9{x; //i) is strictly increasing for 
t G [0,to] and strictly decreasing for t G [to, gl- latter case 9{x; ^i) crosses the value | at 

some do G (0,to) and hence > on [0, do) and n' < on {do, |). This proves the monotonicity 
alternative (al), (a2), and it remains to discuss the number of dislocation eigenvalues. 

We may suppose that the second Neumann eigenfunction is strictly decreasing on [0, d] with its 
unique zero at ^. Thus 9{t; V2) ranges within [|, ^] with 6'(|; V2) = tt, 9'{^; V2) = 1. Therefore, 9 



increases near d/2 and taking into account Lemma 3.10[ a), we find that 9 must be strictly increasing 



on [0,d] and hence R{t;i'2) is strictly decreasing in t on [0,d/2) and on {d/2,d] with 

ii:(0+; z^2) = 0-, f ^-;i/2 j = -00, R('^+;u2j = +00, R{d-; 1^2) = + . (3.21) 
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Case (al): We may suppose that the first Dirichlet eigenfunction is strictly increasing and positive 
^] and even around ^. In this case 6{t; ^i) ranges through [0, |] for t G [ 0, ^] with 6{d/2, jji) = 
[|,7r] for t G [^,d]. Hence 9'{d/2; fii) > 0, and again Lemma 3.10 ^a) implies that 9 



on [0, 

I and through 

is strictly increasing on [0,d] and hence R{t; jJLi) is strictly decreasing in t with 



-R(0+;/xi) = +00, 



; Ail 



0+, 



R{d—; fii) = — oo. 



3.8 



there is no pole At (i.e. 



For t G (0, ^) we have R{t; U2) < < R{t; fii) and hence by Lemma 
St = 0) and there exists a value A G (f2,/^i) with R{t; X) = 0, and this is the only zero. Thus 
we have the uniqueness of the dislocation eigenvalue. For t = |, the zero appears at A = /^i 
which is not inside the gap, i.e., there is no dislocation eigenvalue. Finally, for t G we have 

R{t;i'2) > > R{t; Hi) and hence by Lemma 3.8 there exists a value At G (z^2i/^i)i where R{t; X) 



has a pole. No further poles or zeros can exists, which shows again uniqueness of the dislocation 
eigenvalue. 

Case ( a2): We may suppose that the positive Dirichlet eigenfunction is strictly increasing on [0, do], 
strictly decreasing on [do, 5] and even around i. In this case 6{t; Hi) has the following properties: 



increasing from to ^ for t G [0,do], 
decreasing from 9* to | for t G [to, 



increasing from | to 9* for t G [do, to], 
decreasing from | to 6** for t G [i,d — to], 



increasing from 0* to | for t G [d — toi — do], increasing from | to vr for t G [d — do, d] 
for some to G (do, d/2), which translates into the following behavior of R{t; fii)'- 



decreasing from +00 to for t G [0, do], 
increasing from R^, to for f G [to, f]. 



decreasing from to R^ for t G [do , to] , 
increasing from to i?* for t G [|, d — to] 



decreasing from R* to for t G [d — to, d — do], decreasing from to —00 for t G [d — do, d]. 
If we combine this information with ( |3.21[ ), we conclude: 



(i) For t G (0, do) we have i?(t; 1^2) < < R{t; fj,i) and hence by Lemma 3.8 there exists a value 
A G {i^2,fJ'i) with R(t; X) = 0. No other zero or pole can occur, which shows the uniqueness of the 
dislocation eigenvalue. 

(ii) For t = do the zero has moved to the right-end of the gap, i.e., i?(do; /ii) = 0. 

(iii) Next, we claim that 

R{t; U2) < R{t; ni) < for all t G (do, d/2) (3.22) 
This is obvious for t near do and has to hold by continuity for all t G (do, 5) since equality is 



excluded by Lemma 3.8 Moreover, (3.22) also implies that there cannot be a pole of R{t; X) for 



A G {1^2, l^i) by Lemma 3.8, Thus, R{t; X) increases continuously from R{t; 1^2) to i?(t; ni) as A runs 



through {1^2, f^i) with no zero or pole, i.e., there is no dislocation eigenvalue for t G [do 



' 2^ 



-00, R{d/2; Hi) = and 

R{d/2; X) is strictly increasing for A G [z/2,//i]. In fact, t = ^ leads to a perfectly periodic V{x) due 



(iv) For t = 2 dislocation eigenvalues are excluded since R{d/2; U2+) 



to the symmetry of Vo{x) about x 

'd 



As a result, t = ^ is no dislocation. 



do). For such t we claim that R{t; V2) > R{t; Hi) > 0. Whereas 
positivity is obvious, the ordering is clear for t near |, cf. (3.21), and has to hold by continuity for 



(v) Next we consider t G ( 5 , d 
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all t E {i,d — do) since equality is excluded by Lemma 3.8, Hence, there exists a pole Aj G {1^2, fJ-i 



of R{t; X) and also a zero Aq G {Xt,i^i), which yields exactly two dislocation eigenvalues for t G 
(id-do). 

(vi) For t = d — do, the previous argument still shows the existence of a pole, but the zero has 
moved to the right end of the interval (1/2,^^1) leaving us with only one dislocation eigenvalue. 

(vii) Next, consider t £ (d — do, d). For such t we have R{t; 1^2) > > R{t; /ii) which forces the 
existence of a pole at some value Ai G (2^2, /^i) with no further poles or zeros, i.e., there is exactly 
one dislocation eigenvalue. 

(viii) Finally, t = d is the same as t = and corresponds to no dislocation and hence there are 
no eigenvalues. This completes the verification of the number of dislocation eigenvalues. □ 

Finally, we give a partial answer to the question which of the cases (al), (a2) or (bl), (b2) for a 
given potential Vq actually occur. The condition given in the next theorem is a sufficient condition 
on the potential Vq for (a2), (b2) to occur. 

Theorem 3.13. Suppose Vq satisfies the basic assumptions, i.e., it is even, d-periodic, and con- 
tinuous. 



(i) Assume that Vo is strictly increasing on [0,d/2] and 

'2x 



Vo{x) < V{x) :=/? + (a-/3) 



d 



1 



for allxG [0,d/2], 



where /? := Vo(^) and a G M is arbitrary. If 



[13 - a)d^ > 80(13 - 2V37) ^ 66.75, 



(3.23) 



then only the case (a2) of Theorem 3.12 occurs, i.e., the first Dirichlet-eigenf unction on [0,d] 



is even around ^ but changes its monotonicity at some do G (0, ^). 
(ii) Assume that Vq is strictly decreasing on [0,d/2] and 

Vo{x) < V{x) ■=(5 + {a-p) ^x^ for all x G [0, d/2] , 
where (3 := Vb(0) and a gM. is arbitrary. If 



{p - a)d^ > 80(13 - 2V37) ^ 66.75, 



(3.24) 



then only the case (b2) of Theorem occurs, i.e., the second Neumann-eigenf unction on 
[0,d] is odd around ^ but changes its monotonicity at some do G (0, ^). 



Remark. It will become clear from the proof that (3.23) and (3.24) are not the only conditions 



that lead to the conclusion of the theorem. In fact, by choosing different upper bounds V and a 
different candidate function w{x) in the proof below, one may obtain sufficient conditions which 



are different from (3.23) and (3.24). Since there are manifold ways to derive such conditions, we 



decided to give only the simplest one. Nevertheless, (3.23) and (3.24) are already sufficient to cover 
example potentials such as Vo{x) = sin^(7rx/10) and Vo{x) = cos^(7r2;/10), respectively. 
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Proof. Suppose Vq is increasing on [0, |]. Let fii be the first Dirichlet eigenvalue on [0,d] with 
corresponding positive eigenfunction u. Then ni = kdn, where kdn denotes the first eigenvalue 
on [0, 2] with Dirichlet boundary condition at and Neumann boundary condition at |. Let 6 be 
the Priifer-angle for u normalized by 9{0) = 0, which implies 0{^) = ^. The non-monotonicity of u 
can be shown by proving that e'{^) < 0, i.e., ^ > I in a left-neighborhood of |. By the differential 
equation for 6 we obtain 

e'il) = 1 + (I) - 1) sin2 e =fii-Vo ii) . 

Using the variational characterization of ^1 = kdn, it suffices to find one function w e H^{0,^) 
with w{0) = such that 



/o^ w'^ + Vo{x)w^ dx 



< Vo ii) = 13. (3.25) 



Jq w"^ dx 

Using the upper bound ^0(2;) < V{x) and the quadratic candidate function w{x) = x(2c — x) with 
c G M to be determined, condition (3.25) amounts to 

/ 2 2 — , ^ 2 1^2 

w + V[x)w dx — (3 w dx 
Jo Jo 

= (56(120 -7)c^- 14(240 -7)c + (560 - 7)(i^) <0, where 7 :=(/?- ald^. (3.26) 

3360 

If 7 > 120, then ( |3.26 ) can always be achieved by an appropriate choice of c. If 7 < 120, then the 

d(7-240) 
8(7-120) 

^3(7^ - 2O8O7 + 134400) 



optimal choice for c is c = f[^_?2ni hence (3.26) amounts to 



8 • 3360(120 - 7) ^ °' 



which is fulfilled for 7 G (80(13 - 2\/37), 120) « (66.75, 120). Altogether, the statement (i) of the 
theorem holds true for 7 = (/? — a)d'^ > 80(13 — 2\/37). This concludes the proof of statement (i). 
Part (ii) can be obtained from part (i) by reflecting the interval [0, |]. □ 



Numerical Results We present results of numerical computations of the point spectrum of L 
with the dislocation interface ( |1.3[ ) with s = —t and Vq = sin^(7rx/10) as well as Vq = cos^(-7rx/10). 

As one can see in Figure [7] bottom, the number of eigenvalues in the semi-infinite gap agrees 
Regarding eigenvalues in the firs t finit e gap Gi = {s2, S3) 

(with P 



with Theorem 3.11 

sin^(7rx/10) satisfies the conditions of Theorem 



Vr 



3.13 



note first that since 
1 and, for instance, a = 0.3), 
we know that the first Dirichlet eigenfunction changes monotonicity at some do G (0,d/2) = (0,5). 
The case (a2) of Theorem 3.12, therefore, applies. We obtain numerically do ~ 2.16, see Fig. [Tjtop. 
The number of eigenvalues in the gap Gi agrees with the theory at each t G (0,d), see Figure [7] 
bottom. Eigenvalues in the gaps G2 and G3 are also plotted; note that for these our analysis 



provides no explanation other than the statement of Corollary 3.9 



Figure |8] shows the eigenfunctions corresponding to the 9 labeled eigenvalues in Figure [7j 



The results for Vo{x 
Theorem 



cos^(7rx/10), as an example of a potential that falls in the case (b) of 



3.12, are, in fact, contained in the lower part of Figure [t] because cos^(7r(x — t)/10) 



sin (7r(x — (t + 5))/10). As cos (vrx/lO) satisfies the conditions of Theorem 3.23 (with [3 = 1 and 
for instance, a = 0.3), we know that the alternative (b2) has to apply. 
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02468 10 02468 10 02468 10 



= 10 



t d/2 



do 2 




i 
/ 



0.5 s ^ 



Figure 7: top: the first three band edge Bloch functions of Lq with Vq{x) = sin^(7rx/10); bottom: 
point spectrum of L for (1.3) with s = —t and Vq = sin^(7rx/10) for t G [0,d). The spectral bands 



of L are shaded. Eigenfunctions for the labeled points are plotted in Figure M 
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Figure 8: Eigenfunctions corresponding to the 9 labeled eigenvalues in Figure [Tj 
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3.2.2 One-sided Dislocations 



As the second representative example of the dislocation problem (1.1) with (1.3) we choose s 
0,t £ (0,d), which is equivalent to the system 



^0 := -d2^t ^ Vo{x)'^^ = Xi)^ for x < 0, 
L*.^* := -dlil)^ + VqIx + t)ip^ = A^* for x > 



(3.27) 



coupled by the the C^-matching conditions (3.15). Localized eigenfunctions ip^, once again, exist 
only for A ^ cr(Lo) and have the form 

i^\x; A) = X{x<o}i^-{x; ^) + X{x>o}V'+(2;; A), 

where 'tp^{x; A) and V'+C^! ^) ^-re those Bloch functions of Lq and Lt = —d'^ + Vo{x + t), which decay 
on R~ and M+, respectively. The matching condition ( |3.15 ) now becomes 



Rlit; A) = i?0(0;A), 



where R^{0; A) is the same as i?_(A) defined in (3.4). 

Because R^{0; A) is decreasing and continuous in each gap (Lemma 
of A) in Section 



3.2.1 



3.1) and given the analysis 
determining intersections of R^{t] A) and RzJO; A) in Go and Gi is now 



st r aight forward . 

Lemma 3.14. For s = the number of dislocation eigenvalues in any gap Gn,n > 0, is 0, 1 or 2. 

Proof. R\{t; A) is strictly increasing and continuous in A on each continuity segment and its conti- 
nuity can be broken only at one point (pole) in G„, see Lemma 



3.8 



As R^ (0; A) is continuous and 



decreasing throughout G„, only up to 2 intersections of R^{t; A) and R_{0; A) can occur. 



□ 



Theorem 3.15. Suppose Vq satisfies the basic assumptions, i.e., it is continuous, even and d- 



periodic, and let s = in (1.3), and consider the semi-infinite gap Gq = (— oo, si). 



(a) If Vq is strictly increasing on [0,d/2], then there is no/exactly one dislocation eigenvalue in 
Go for t £ [0,d/2] / {d/2,d) respectively. 

(b) If Vq is strictly decreasing on [0,d/2], then there is exactly one/no dislocation eigenvalue in 
Go for t £ (0,d/2) / [d/2,d] respectively. 

Proof. We, once again, present the proof only of (a) as (b) follows by shifting the potential in x 



(or t) by d/2. As explained in the proof of Theorem 3.11, si is a Neumann eigenvalue and the 
corresponding eigenfunction can be taken positive on [0,d/2] with n' < on (0,d/2) and with a 
point of even symmetry at x = d/2. 

the function (0; A) decreases continuously from oo at A ^ — oo 

It follows 



By Lemmas 



3.1 



3.2 



and 



3.3 



3.11 



to at A = si. The behavior of R^{t; X) is explained in the proof of Theorem 
that R^{0; X) and R^{t; X) intersect in Gq exactly once for t E {d/2,d) and do not intersect for 
t£[0,d/2]. □ 



Theorem 3.16. Suppose Vq is an even, d-periodic C"^ -function, let s = in (1.3), and consider 
the first finite gap Gi = (s2, S3). 
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(a) If Vq is strictly increasing on [0, 2]; and hence Gi = (z/2,/Ui), then there is exactly one 
dislocation-eigenvalue in Gi for all t G (0,d). 

(b) If Vq is strictly decreasing on [0, and hence Gi = (/Ui,z^2), then we have the following 
alternative for the second Neumann-eigenf unction u: 

(bl) u is strictly monotone on [0, Then there is exactly one dislocation- eigenvalue in Gi 
for all t G (0,d). 

(b2) u changes monotonicity on [0, exactly once at the extremal point do G (0, |). Then 
the number of dislocation-eigenvalues in Gi is as follows: 



dislocation parameter 


t G (0,do) 


t G [do, d - do) 


t G [d-do,d) 


number of eigenvalues 


2 


1 






Proof. Case (a): As explained in the proof of part (a) of Theorem 3.12 for t G (0,d/2] we have 



R^{t;v2+) < and R^{t; X) continuous and increasing in A G Gi. Therefore, R^{t; X) intersects 
R^{0; A) exactly once, as R^{0; A) decreases continuously from at A = 1^2+ to —00 at A = fJ-i—, 
see Lemmas 13.11 13.21 



Next, as the proof of Theorem 



3.12 



(a) shows, for t G {d/2, d) the function R^{t; A) has a pole at 
some Xt G Gi and increases continuously on the interval (^2, Xt) with R%{t; 1^2) > 0, R^{t; Aj— ) = 00 
and on the interval [Xt, ^1) with Ai+) = —00. The functions A) and R^{0; A), therefore, 

intersect exactly once on A G {Xt,fJ.i) and they do not intersect on A G {^21 A*). 

Case (b): In the case of Vq strictly decreasing on [0,(i/2] the function R^{0;X) is continuous and 
strictly decreasing from 00 at A = 1/2+ to at A = /ii, see Lemm as |3.H |3.2[ We obtain below the 

by the shift of d/2 in t 



3.12 



(a) 



behavior of A) from that of R{t] A) in the proof of Theorem 

and switching of the roles of and 1^2. 

Case (bl): For t G (0,d/2) we have Rlit; fii) > > Rl{t;i^2) and Rlit; X) has one pole in 
A within Gi. R^{0; X) thus intersects R^{t; X) exactly once on Gi. For t G [d/2,d) the function 
R+it] A) is continuous on Gi and R\{t; fii) < < R^{t; 1^2— )■ Exactly one intersection of R\{t; A) 
and R^{0; A) thus exists. 



Case (b2): For t G (0,(io 



in (v) in the proof of Theorem 3.12 (a2). 
proof of Theorem 



3.12 



the function R\(t; X) behaves in A like R{t] A) on t G {d/2, d — do) 
Note that do here corresponds to d/2 — do in the 



Namely, we get Rl{t; m) > Rl{t;u2) > and a pole of i?°(t;A) at 
some A( G Gi. Two intersections of R^{t; X) and i?^(0; A) thus exist. For t G [do,d — do) the 
behavior of the eigenfunction and hence of R^{t; X) is qualitatively the same as in (bl) of this 
proof and precisely one eigenvalue thus appears in Gi. Finally, for t G [d — do,d) the function 

(a2). 



3.12 



R^{t] X) behaves in A like R{t; X) on t G [do, d/2) in (iii) in the proof of Theorem 
Therefore, R^{t; ni) < i?5]_(t;z^2) ^ and R^(t; X) is continuous throughout Gi. No intersections 
of Rl{t;X) and R^_{0;X) thus occur. Finally, t = d corresponds to no dislocation resulting in a 
purely continuous spectrum of L. □ 



Numerical Results 

(1.3) with s = and Vq 



and 



Results of numerical eigenvalue computations with the dislocation interface 
sin^(7r3;/10) are displayed in Figure k)l They agree with Theorems 3.15 



3.16 Figure 10 shows the eigenfunctions corresponding to the 6 labeled eigenvalues in Figure|9| 
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As expected, they lack symmetry in contrast with the eigenfunctions of the symmetric dislocation 
in Figure [8} 




1 (t = 8.2, A » 0.22) 



2 (t = 9, A«0.66) 



3 (t = 2.5, A « 0.49) 




4 (t = 8, A»0.96) 



-20 20 

X 

5 (t= 1.6, A« 0.98) 




•20 20 

X 

6 (t = 5, A«1.42) 



0.02 
0.01 



-0.01 
-0.02 



-50 



-400 -200 200 400 



Figure 10: Eigenfunctions corresponding to the 6 labeled eigenvalues in Figure [9] 

The results for Vb(x) = cos^(7rx/10), as an example of a potential that falls in the case (b) 
of Theorem |3.16[ appear in Figures [TT] and 12 As we know from Section 3.2.1 the potential 
cos^(7rx/10) falls into the case (b2) and the second Neumann eigenfunction thus changes mono- 
tonicity on (0,(i/2), see Figure [TT] top. Agreement of the numerics with Theorems 3.15 and 3.16 is, 
once again, observed. 
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supported by the Alexander von Humboldt Research Fellowship. 
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Figure 11: top: the first three band edge Bloch functions of Lq with Vo{x) = cos^(7rx/10); bottom: 
point spectrum of L for (1.3) with s = and Vq = cos^(7rx/10) for t G [0,d). The spectral bands 
of L are shaded. Eigenfunctions for the labeled points are plotted in Figure [12} 




Figure 12: Eigenfunctions corresponding to the 6 labeled eigenvalues in Figure 11 
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